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ABSTRACT 

A number of Kepler planet pairs lie just wide of first-order mean motion resonances 
(MMRs). Tides have been frequently proposed to explain these pileups, but it is still 
an ongoing discussion. We contribute to this discussion by calculating an optimistic 
theoretical estimate on the minimum initial eccentricity required by Kepler planets 
to explain the current observed spacing, and compliment these calculations with N- 
body simulations. In particular, we investigate 27 Kepler systems having planets within 
6% of the 2:1 MMR, and find that the initial eccentricities required to explain the 
observed spacings are unreasonable from simple dynamical arguments. Furthermore, 
our numerical simulations reveal resonant tugging, an effect which conspires against 
the migration of resonant planets away from the 2:1 MMR, requiring even higher 
initial eccentricities in order to explain the current Kepler distribution. Overall, we find 
that tides alone cannot explain planets close to 2:1 MMR, and additional mechanisms 
are required to explain these systems. 


1 INTRODUCTION 


The NASA Kepler mission has been immensely successful for 
detecting planets outside of our solar system. To date, it has 
discovered over 4,500 exoplanet candidates along with 466 
multi-planet systems ( [Akeson et al.||2013] [Rowe et al.|2014| ) . 
A number of these systems are just wide of a mean motion res¬ 
onance (MMR), which occurs when the period of one planet 
is an integer ratio of another. In particular, statistical excesses 
in the period distribution of Kepler planets have been detected 
just wide of the 3 :2 and 2:1 MMR (iLissauer et al.[|2011[|Fab-| 


|rycky et al.|2014| [Steffen & Hwang|2015 ). It is believed that 

in the past these planetary systems migrated into resonance 
via convergent migration ( [Lee & Peale||2002| ), and a num¬ 
ber of dissipative mechanisms have been proposed to slowly 
bring these planets out of MMR. The most popular dissipative 
mechanisms to explain the observed near-resonant systems 
are tidal ( Lithwick & Wu|[20T2l [Batygin & Morbidelli||2Q13[ 


Delisle et al.[[2Q14 ), protoplanetary ( Rein[[2012[ [Baruteau &[ 

Papaloizou[[2Q13[ Goldreich & Schlichting[|2014P, an d plan- 

etesimal ( Moore et al.|2013 Chatterjee & Ford|2()15| ). 


In this work we focus exclusively on tidal dissipation, for 
which there is no clear consensus on whether this mechanism 
alone can successfully explain the excess of near-resonant 



Figure 1. Kepler systems close to 2:1 MMR. A statistical excess is 
present just wide of the 2:1 MMR, and appears to decline beyond 6% 
of the resonance, as marked by a red dotted line. 


pairs near first order MMRs. Several authors ([Lithwick & Wu 

L 

2012[ Batygin & Morbidelli[[2013[) have argued for (Delisle 

: used this relationship to show that most near-resonant planet 

et al.[2014) tidal dissipation, whereas others ([Lee et al. 

2013) 

pairs cannot be explained via this mechanism, the few ex¬ 
ceptions being small rocky planets for which tidal dissipation 

have argued against it. In particular, [Lithwick & Wu ( 

2012) 


first introduced the mechanism of resonant repulsion, and 
showed that in the limit of low eccentricities for near-resonant 


planets, the space between planets grows as Batygin & 
|MorbidelII| ( [2Q13| ) confirmed this result. [L^ et al.[ ( |^13| ) then 


a high eccentricity mechanism by which planets may still be 
able to evolve to their current positions via tides alone. 

We contribute to this debate by developing optimistic 
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2 Silburt & Rein 


theoretical estimates for the evolution of planets away from 
resonance, and compare these estimates to N-body simula¬ 
tions. We then make a statement about the likelihood of tides 
to explain near-resonant pairs. We focus on Kepler systems 
within 6% of the 2:1 MMR, which appears to be the natural 
cutoff for this statistical excess, as shown by a dotted line in 
Figure Using this sample we will calculate the minimum 
required initial eccentricity to explain their current positions, 
given that they started in MMR and evolved under the influ¬ 
ence of tides alone. In addition, we present numerical findings 
of resonant tugging, an effect which prevents the evolution of 
planets away from MMR when eccentricity is high, making it 
significantly more difficult to achieve the observed spacings 
seen today. Resonant tugging does not appear to have been 
extensively studied/accounted for in this field due to the fact 
that most analysis of planets in resonance have worked in the 
e <C 1 limit, and we find that resonant tugging exclusively 
affects planets in MMR with moderate to high e. 

Our paper is organized as follows - in Sectionj^we out¬ 
line our theoretical and numerical framework, in Section 
we present the main findings of this paper, in Section]^ we 
present a discussion and conclude in Section 


2 METHODS 


2.1 Theory 


The equations widely used to describe the evolution of plan¬ 
ets under the influence of tides are (e.g. |Barnes et al.|2008t 
iLithwick & Wu|2Q12t[Lee et al.|2013l : 


9 k 1 

IgM^ frp\^ 

V a3 [a) " 

(1) 


Q rUp \ 


(2) 


where a is the semi-major axis, e is the eccentricity, k is the 
planet’s Love number, Q is the planet’s tidal quality factor 
( [Goldreich & Soter||1966| ), rup and Vp are the planet’s mass 
and radius, respectively, M is the stellar mass, and G is the 
gravitational constant. From Eqs.[^and[^it follows that e and 
a are related by: 


- = 2ee (3) 

a 

Which arises from conservation of orbital angular mo¬ 
mentum. We are interested in finding a relationship between 
the total migration of the system and its state variables (such 
as a and e), which can be obtained by integrating Eq.[^ 

r^^ = 2rede 

Jaf Jcf 


Of = Ui exp(-ei + ej), (4) 

where subscripts i and / refer to the initial and final states, 
respectively. 

Eq .[^is surprisingly simple - if we know the initial states 
Ci and Ui of a planetary body, as well as the final eccentricity 
e/, we can predict its final position. Eq. [^is independent of 


tidal response parameters {k, Q), the length of time consid¬ 
ered, the mass and radii of the star/planet, etc. These addi¬ 
tional factors affect the timescale by which the body arrives 
at its final position, a/, but not the final position of the body 
itself. 

Let us measure the spacing of a planet pair close to a 
j : j + 1 MMR by defining 


_ Pout J 1 

Pin j 


(5) 


where P is the orbital period, and the subscripts in and 
out refer to the inner and outer planet, respectively. For the 
2:1 MMR, j = 1. This definition of A is the same as Lee et al.| 


( [2Q13| ), and twice the value used by |Lithwick & Wu 


_ __ 

Using the fact that P oc and substituting Eq7|4[Into 

Eq. we recast A in terms of ai and for the near resonant 
pair (setting j = 1): 

3/2 


A = 


Uoutjf 

ain,f 


- 2 


( 6 ) 


_ ( ^out,i exp( ContjZ ^out,f^ 

\ G^^n,zexp(-e2^ . 

We make the assumption that eout,i ~ (see Section 

for a discussion). After simplifying, we get: 

Uoutji ^in,f^ 

ain,i 

Thus, Eq.[7| relates the final spacing of the system. A, to the 
initial conditions of the system (a^ and ei) and the final ec¬ 
centricity of the inner planet, Cinj. Eq. [^assumes that the 
angular momentum of each individual planet is conserved, 
which in general is not true for multi-planet systems (only 
total angular momentum is). Other works ( [Lithwick & Wu| 
|2012[ [Batygin & Morbidelli||2013] ) have accounted for this 
fact, resulting in more accurate equations for the evolution 
of a planet pair. We now aim to compare Eq. to numerical 
simulations, and estimate its accuracy for Kepler systems with 
A < 0.06 of the 2:1 MMR. We first outline our experimental 
setup. 





2.2 Experimental Setup 


Numerical simulations are carried out using the Wisdom & 
Holman integration scheme ( Wisdom & Holman||1991| ), im¬ 
plemented via REBOUND ( [Rein & Liu|2012 ). Our sample con¬ 
sists of Kepler systems near the 2:1 MMR with A < 0.06, 
which we call near-resonant pairs. Our choice of A < 0.06 
stems from a natural cutoff where the excess of near-MMR 
ends, as shown in Figure We exclude near-resonant pairs 
interior to the 2:1 MMR, since tidal forces appear only to in¬ 
crease planet separation with time (e.g. [Lithwick & Wu[2012| ). 
In addition, we also exclude near-resonant pairs that are com¬ 
plex. By complex, we mean dynamically involved in an ad¬ 
ditional (near) resonance (e.g. Laplace resonance), and/or 
containing an additional planet orbiting between the near¬ 
resonant pair. Our exclusion of complex resonant systems de¬ 
creases the number of Kepler systems by 6. Lastly, we also re¬ 
move the Kepler-11 system, which does not contain a complex 
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To speed up simulation time, we increase /c/Q by a fac¬ 
tor of Ak/Q (or alternatively this could also be interpreted as 
increasing tidal strength). This tactic has been used by other 
scientists (e.g. |Delisle et aL|2014| ), and is valid as long as Te is 
much longer than the planet’s eccentricity libration time. We 
simulate our Kepler sample for 50 Myr, and use A^/q — 200, 
giving a total simulation time of T = 10 Gyr. 

To begin our simulations, we place each planet at a dis¬ 
tance of lASttobs, where aobs is the observed semi-major axis 
value from the Kepler catalog, and assign a = 0.01. We then 
migrate each planet (in a type-I fashion) back to its original 
starting position aobs except for of the outer near-resonant 
planet, which instead migrates a distance of aobs + A, forcing 
the near-resonant pair into a 2:1 MMR. 

Each planet migrates for time tmig at rate a = anp!^^^ jC, 
where n is the mean motion of the inner planet, p is the 
planet/star mass ratio, and C is a constant. Lower values 
of C cause the outer planet to encounter the MMR sooner, 
allowing time for both planets to migrate together in reso¬ 


nance, which increases eccentricity to a desired value (Lee 


& Peale |2002|). F or the restricted 3-body problem |Goldre- 


ich & Schlichtingl ( [2014| ) guarantee capture into resonance 
if Cout > 3.75, Cin < {pout/jY^^ but we use a more conser¬ 
vative value of Cout = 6 as well as perform numerical tests to 
ensure that overstabilities do not occur on Myr timescales. 

Defining K = , we use a default K — 100 when 

migrating planets into resonance but also experiment with 
K — 10. iT (along with affects the resulting equi¬ 

librium eccentricity, but does not affect tidal evolution, and 
thus does not affect our main conclusions. At this point, ini¬ 
tial eccentricities of our simulated Kepler systems range from 
0.05 < Ci < 0.25, depending on the value of C, K, etc. 

After time tmig, migration is quenched over a timescale 
of tmig IS by letting Ta ^ oo and Te oo, where Ta = —a/a 
and Te = —e/e. It is at this point that tides are turned on, and 
the system evolves under the influence of Eq.[^and|^for the 
remainder of the simulation. 


tions and the observed Kepler spacing. 

3 RESULTS 

3.1 Theory vs. Numerics 

We compare our theoretical predictions of tidal evolution, 
^th (Eq.|7j), to our numerical simulations, ^Yium (Eq. The 
difference, Anum — Ath, is displayed as a solid line in Fig¬ 
ure and expressed as a cumulative distribution (GDF). For 
all but 2 systems we see that Anum — Ath < 0, indicating that 
our theoretical predictions Ath consistently over-predict our 
numerical results, Anum- These 2 exceptions, Kepler-32 and 
Kepler-221, are expected due to the fact that T/re ~ 1000, 
which allows extensive A oc t^^^ resonant repulsion growth 
dLithwick & Wu|201^ in time T which is not accounted for in 
Eq .[t] Otherwise, we find that Eq.[7|consistently overestimates 
the true evolution of the system. 

The reason for this Anum - Ath < 0 trend is due to 
resonant tugging, an effect present in the numerics but not 
captured by our theoretical predictions. Resonant tugging acts 
to keep planets closer together than theory would predict. We 
explore resonant tugging in the next section. 

For the same simulations we also plot Anum — Aobs as a 
dashed line in Figure which is the difference between our 
numerical results, Anum, and the observed spacing of Kepler 
planets seen today, Aobs (Eq.[^. As is clearly shown, Anum — 
Aobs < 0 for all but two systems, suggesting that tides alone 
cannot explain the observed spacing of Kepler planets. These 
two exceptions are again, Kepler-32 and Kepler-221, and are 
exceptions for the same reason as above. Since Aobs = 0.038 
and 0.035 for Kepler-32 and Kepler-221, respectively, using 
a lower A cutoff for our Kepler sample (e.g. Aobs < 0.03) 
would not have changed our result that tides cannot explain 
near-resonant pairs. 

Although highly suggestive, this result does not conclu¬ 
sively disprove tides as the primary evolving mechanism, 
since higher a could cause more migration in time T (see 
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4 Silburt & Rein 


Eq.il, and the median eccentricity for our simulations is 
ein,i = 0.14. Higher values of ein,i are possible. Instead of 
numerically exploring every possible ein,i value, however, we 
instead reverse the argument in Section [3^ and calculate the 
minimum ein,i required to explain Aobs- First, however, we 
explore resonant tugging. 

3.2 Resonant Tugging 

Resonant tugging affects planets in MMR subjected to en¬ 
ergy dissipation (e.g. tidal), with moderate to high eccen¬ 
tricity. When these conditions are met and the inner planet 
migrates inwards (trying to leave the resonance) it tugs the 
outer planet inwards along with it, transferring dissipative 
forces from the inner to the outer planet. The result is that 
the inner planet migrates less than expected, the outer planet 
migrates more than expected, and the planets are closer to¬ 
gether than theory would have predicted. 

Figure illustrates resonant tugging for a pair of m = 
10“^M planets in MMR. For the black curve ein,i = 0.125, 
for the grey curve ein,i = 0.018, and otherwise the initial 
conditions of each test case are the same. In both cases we 
allow only the inner planet to evolve under the influence of 
tides. The top and bottom panels show the period evolution 
of the outer and inner planets, respectively, and the dotted 
curve in the bottom panel shows evolution of the inner planet 
in the absence of the outer planet (also ein,i = 0.125). 

Resonant tugging is exhibited in the first 0.5 Gyr of evo¬ 
lution for the black curve, showing how tidal forces affecting 
the inner planet also affect the outer planet by dragging it 
inwards too. Comparing the solid and dotted black curves in 
the bottom panel of Figure we see that in the presence of 
the outer planet, the inner planet migrates much less than ex¬ 
pected due to resonant tugging. Since the outer planet has 
also migrated inwards more than expected, the result is that 
Abiack has grown very little over the first 0.5 Gyr of evolution 
{Abiack = 0.004 after 0.5 Gyr). 

Within the framework of resonant tugging, the outer 
planet can be thought of as a massive anchor - as 
moutimin oo the inner planet has an increasingly diffi¬ 
cult time migrating both bodies inwards, leading to a pair of 
(relatively) stationary planets. Conversely, as moutimin 0 
the inner planet has an easier time migrating both bodies in¬ 
wards, and the trajectory of the inner planet will approach its 
single-planet trajectory (i.e. the dotted black line in Figurej^. 

We now briefly contrast resonant tugging from resonant 
repulsion (first described by |Lithwick & Wu]|2Q12| ). Qualita¬ 
tively, the two main differences between resonant tugging and 
resonant repulsion are: 

(i) Resonant tugging exclusively affects planets in reso¬ 
nance (resonant angle (s) librating) with moderate to high ec¬ 
centricity, while resonant repulsion affects planets both in and 
close to resonance, and is most noticeable in the e ^ 1 limit. 

(ii) Resonant tugging decreases Uin and decreases Uout, 
while resonant repulsion decreases ain and increases Uout- 

These two differences are illustrated in Figurej^ The grey 
curve, which has low initial eccentricity (ein,z = 0.018) ex¬ 
hibits pure resonant repulsion - a decrease in ain and an in¬ 
crease in Uout, while the black curve first exhibits a shorter 



time (Gyr) 

Figure 3. Two test cases illustrating resonant tugging and repulsion. 
The top and bottom panel shows the period evolution of the inner 
and outer planet, respectively. For the black curve Cin^i = 0.125, 
while for the grey curve ein^i =0.018. The dotted black curve shows 
the numerical trajectory of the inner planet {cin^i = 0.125) in the 
absence of the outer planet. 

period of resonant tugging followed by resonant repulsion. 
The transition from resonant tugging to resonant repulsion 
for the black curve occurs after 0.5 Gyr, when the eccentricity 
of the inner planet has dropped to a low (e^n = 0.035) value. 

In Figure [U after ~ 1 Gyr, Auack < Agrey, but Agrey = 
Abiack, showing how resonant tugging can permanently stunt 
the growth of A. We also see for the black curve that aoutj < 
aout,i, and it is only after many Xe timescales that the outer 
planet can recover (or exceed) its initial position via resonant 
repulsion. Furthermore, since to first order (when T/xe 
1000) A oc Ta,in oc One would naively expect the black 
curve to experience more A growth in time T, yet we ac¬ 

tually find Abiack < Agrey, illustrating just how significant 
resonant tugging can be. This means that the high eccentric¬ 
ity tidal mechanism suggested by [Delisle et al.| ( [2014| ) does 
not work for planets in resonance, since (contrary to expec¬ 
tation) high eccentricity actually stunts the growth of A, not 
accelerates it. 

The results shown in Fig. are due to resonant tugging, 
and is supported by the fact that for every simulated Kepler 
system we find aoutjluout^i < 1. Since resonant repulsion 
can only increase Oout (as shown in Fig. , OoutJ / Oout^i ^ 1 
can only be due to resonant tugging since inward tidal mi¬ 
gration is a negligible contribution for the outer planet. The 
analytic aspects of resonant tugging will be studied in more 
depth in future works. 


3.3 Minimum Eccentricity to Explain Aobs 

Since we found in Section |3.1| that Ath consistently over¬ 
predicts the amount of tidal migration (due to resonant tug¬ 
ging), we can use it as an upper limit predictor of tidal evo¬ 
lution, assuming that near-resonant pairs started in 2:1 MMR 
and evolved to their present locations. Starting from Eq. 
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Figure 4. Three CDFs showing the theoretical minimum eccentricity 
required by the inner planet in order to achieve the observed A spac¬ 
ing seen by Kepler planets today The solid, dashed, dotted lines repre¬ 
sent T = 1, 5,10 Gyr tracks, respectively, while the dash-dotted line 
represents T —oo. In all calculations we assume the outer planet 
remains stationary, i.e. (ai/af)out = 1- The red shaded region marks 
the unphysical region where the eccentricity is larger than unity The 
blue region marks the region where most systems undergo a dynam¬ 
ical instability 


and using the same logic as Section [3J] we calculate the min¬ 
imum eccentricity required by the inner planet to achieve the 
observed spacing seen today, Aobs, after T years (see Ap¬ 
pendix]^ for detailed calculations): 


^in,i 


\ 


1 - exp(-2T/Te,zn) 


( 8 ) 


Thus, for a pair of planets starting in 2:1 MMR, if we know the 
observed spacing today, Aohs, the number of Te,in damping 
timescales in time T, and estimate the amount of migration 
done by the outer planet in time T, {ai/af)out, we can cal¬ 
culate the minimum initial eccentricity that the inner planet 
must have in order to arrive at the current spacing, ein,i- Since 
we assume a starting position of exact 2:1 commensurability. 
Equation is only tailored for Aobs > 0. Tides cannot de¬ 
crease planet spacing over time. 

Figure[^shows CDFs of Eq. [^applied to our Kepler sample 
for T = 1 (solid), 5 (dashed), and 10 (dotted) Gyr with the 
outer planet remaining stationary, i.e. {aijaf^out — 1 (see 
Section]^ for a discussion). In addition we plot a T ^ oo 
curve as a dash-dotted line, which the other curves converge 
to. The shaded red region marks where Cin.i ^ 1.0, while 
the blue marks an unstable region where eccentricities are 
unlikely to exist. We construct the blue region by numerically 
finding the maximum eccentricity allowed before > 50% of 
our Kepler systems go unstable within 2 Myr (see Appendix]^ 
for details). We find this maximum eccentricity to be 0.3 
(the left boundary of the blue shaded region in Figure]^. 

The three T < oo curves in Figure]^ lie largely in the red 
and blue shaded regions, indicating that the required eccen¬ 
tricities to explain Aobs are unreasonable. We thus conclude 


that most Kepler systems cannot be explained by tides alone. 
Even for the T = 10 Gyr curve, an optimistic estimate for the 
age of many Kepler systems, about 35% of systems still cannot 
be explained due to tides alone. Glearly, another mechanism 
is needed to explain the near-resonant 2:1 MMR pairs. 


4 DISCUSSION 


A number of assumptions have been made in constructing Fig¬ 
ure]^ which only strengthen our conclusion that planets close 
to the 2:1 MMR cannot be explained due to tides alone. 

First and foremost, from Sections [3.1| and [3^ we have 
shown that resonant tugging causes our theoretical predic¬ 
tions of A to be overestimates, and thus Eq. [^underestimates 
the minimum ein,i required to explain the current observed 
spacing. This discrepancy becomes most pronounced when 
moutimin > 1, which is the case for many Kepler systems 
in our sample. In particular we showed that resonant tugging 
can stunt t he ev olution of planets away from MMR by roughly 
(see Sec jSjj ), and this stunted evolution is not accounted 
for in Figure]^ 

Second, we have assumed optimistic k/Q values which 
allow for more migration in time T. Although some planets 
could have such generous values, it is unlikely that all of them 
do. Since a oc {k/Q)e^, as k/Q decreases ein,i must increase 
in order for the planets to achieve the same observed A in 
time T. 

Third, our estimate of Cmax = 0.3 (blue-shaded region 
in Figure]^ is very likely an overestimate. For example, jPu &| 
|Wu| (|2015| ) found that the mean eccentricity of high-multiple 
Kepler planets must not exceed Cmean = 0.02 to guarantee 
long-term dynamical stability. Also, from simple orbit crossing 
arguments when e > 0.23 the perihelion of the outer planet 
crosses the aphelion of the inner planet for a 2:1 MMR, and 
long-term stability can no longer be guaranteed. Even if it 
were possible for Kepler systems to remain stable with high 
(> 0.3) eccentricities, it is unclear what kind of mechanism 
could consistently generate them for our Kepler sample. 

There are other assumptions we made throughout this 
paper which we now justify. For the results presented in Fig.]^ 
we assumed that the outer planet remains stationary, i.e. 
{ai/af)out = 1.00. This is a reasonable assumption since, re¬ 
ferring to our T = 10 Gyr numerical simulations as a bench¬ 
mark, the median value of {ai/af)out = 1.001 1. Increasing 

the value of {ai/af)out shifts our GDFs in Fig. further into 
the blue/red instability region. 

In Eq. [^we assumed that eout,i ~ ^ontj which essen¬ 
tially states that the the initial and final positions of the outer 
planet are the same. As stated in the previous paragraph, 
since we found numerically from our simulations that the me¬ 
dian {ai/af)out ~ 1, this assumption is reasonable. 

In constructing FigJ^we have used Eq.|^ which is Eq. 89 
from |Delisle et al.| ( [2014 ), who argued that for moderate to 
high eccentricities iein,i ^ 0.15), many of the near-resonant 
pairs could in fact be explained by tides. There are a number 
of differences between our analysis and theirs however. First 
(and most importantly), their estimates are based on theoret¬ 
ical predictions (i.e. resonant tugging unaccounted for), and 
we have shown in Section [3^ that the growth of A is signif- 
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icantly stunted when resonant tugging is accounted for, es¬ 
pecially when rrioutlmin > 1. Second, their analysis assumes 
that T ^ oo, while we restrict to T = 1,5 and 10 Gyr. Looking 
at Fig we see that the T ^ oc curve tells a very different 
story than the T = 1,5 and 10 Gyr curves, and the conclu¬ 
sion of whether or not tides can explain Aobs is certainly time 
dependent. Lastly, [Pelisle et al^ assumes A = 0.03 for all sys¬ 
tems, while we use the system specific A values. 

As a consistency check for our results, we perform the 
same set of experiments (including our tests of resonant tug¬ 
ging) using a different version of tides, implementing them in 
terms of forces (as opposed to orbital elements like in Eq. 
and|^ according to |Papaloizou & Larwood] pOOQ| ) : 

(9) 


where a damp is the damping acceleration, v is the velocity, r is 
the position, r is the scalar position and Xe = — e/e as before. 
Whenever a planet receives a “kick” in the WH integration 
scheme, an additional kick of adamp is supplied to account for 
tides. We find our overall conclusions unaffected using this 
implementation of tides. 

We have omitted complex resonances from our analysis 
because their behaviour is much more unpredictable. For sim¬ 
ple resonances we found that Ath > Anum, but some com¬ 
plex resonances violated this relationship. The reasons are 
currently unknown, but will be more thoroughly investigated 
in future works. 

Lastly, it should be mentioned that spin tides were omit¬ 
ted from this analysis, which arise when the spin rate of the 
host star Q* is different than the mean motion n of the orbit¬ 
ing planet (responsible for the Moon’s recession from Earth 
over time). Spin rates of Kepler stars are largely unknown, as 
well as the evolution of these spin rates, dQ^/dt. Depending 
on the sign of (Q>, — n), spin tides can induce inward or out¬ 
ward migration. It is thus a non-trivial process to determine 
what the affect of spin tides might be on the evolution of a 
system. The equation governing spin tide migration is ( [Mur-| 
|ray & Dermott|199^ : 


. . . 3/c* rup f R\ 

ap = sign(f2* -np)^^ ( — ] ripttp 
M \ap J 


( 10 ) 


where the subscripts p and * refer to the planet and star, re¬ 
spectively. We can estimate the relative strength of eccentric¬ 
ity tides (Eq. to spin tides (Eq. [To] ) : 

q _^ / gm3 

decc _ Qv\l Q* f ^ \ 

dspin q 'i^p f ^ ^ ^ n V ^ -R ^ 

M \apj '^P^P 


Interestingly enough, the relative strength of spin vs. eccen¬ 
tricity tides is independent of the semi-major axis. Assigning 
typical values from Wu & Murray ( 2003| ) for Q*, (k/Q)^, M 
and R, and assuming a typical ~ 4m^ planet we get: 


r^ecc 

(^spin 


30 


( 11 ) 


Gombining this with the fact that and dQ.^/dt are largely 
unknown for Kepler stars, we felt justified omitting spin tides 
from our analysis. 


5 CONCLUSION 

In conclusion, we have investigated 27 Kepler systems con¬ 
taining 2:1 near-resonant pairs, and find that tides alone can¬ 
not explain their current observed spacing, Aobs- In Figure 
we calculated the minimum theoretical eccentricity required 
by the inner planet to explain Aobs and found that for a large 
number of systems ein,i > 0.3, which from simple dynamical 
arguments is not a reasonable eccentricity for Kepler planets 
to have. Furthermore, our numerical study of resonant tug¬ 
ging reveals that our theoretical predictions of ein,i are opti¬ 
mistic estimates, and in cases where rrioutlmin > 1, signifi¬ 
cantly so. A number of other assumptions made throughout 
the paper contribute to these optimistic estimates. 

As a numerical compliment to our theoretical investiga¬ 
tion, we simulated our Kepler sample for 10 Gyr with a me¬ 
dian eccentricity of ein,i = 0.14, and found only two systems, 
Kepler-32 and Kepler-221, that migrated to Aobs (dotted line 
in Figure H- Glearly another mechanism is required to ex¬ 
plain the excess of Kepler systems exterior to the 2:1 MMR. 
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APPENDIX A: CALCULATION OF MINIMUM 
ECCENTRICITY REQUIRED FOR INNER PLANET GIVEN 
I^OBS AND T YEARS 

Starting from Eq. we now derive Eq. 

- = 2ee 
a 


da ^ oc , ^^^ 

— = 2e^-dt (Al) 

a e 

Since Xe = —e/e is a constant, the only quantity with a time 
dependence on the right hand side is e. Integrating Eq. [i]g ives 
us an expression for e(t): 


de 

e 


9 

2 Q rup 



e 1 

dt = -dt =- dt 


e Te 


Since (again) we assume Xe = —e/e is a constant, the inte¬ 
gration is straightforward to yield: 


-/l) (A2) 

Now that we have eccentricity as a function of time, we 
can plug Eq. |A2| into Eq. |Al| and integrate to get: 



e{t) = Ci exp 


ln(ai/a/) = ef ( 1 — exp(—2 — T) 

'l~e 


We rearrange Eq.l^for a to get: 


\n{ai/af) 


1 - exp(-2T/xe) 


(A3) 


(A4) 


Where all quantities refer to a single planet (e.g. the inner 
planet). To connect Eq.j^to a 2:1 MMR pair, we plug in Eq.[^ 
along with the fact that ain,i — (for 2:1 MMR) to 

get: 


\n(ai/af)i 


1 - exp(-2T/xe,zn) 


S 

fN 

C 
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Figure Bl. The fraction oi Kepler systems in our sample that go un¬ 
stable within 2Myr if each planet is given an initial eccentricity of . 
The error bars derived from Poisson statistics. 


^in,i 




In rr \ ( (^obs+2)^/^ 

22/3 )( f 

1 - exp(-2T/xe,zn) 


^in,i 




1 - exp(-2T/xe,zn) 


which is the result displayed in Eq. 


(A5) 


APPENDIX B: MAXIMUM ECCENTRICITY FROM 
DYNAMICAL SIMULATIONS 

In constructing the blue shaded region in Figure]^ we have 
performed numerical simulations of stability for each Kepler 
system. We simulate each Kepler system for 2 Myr, and assign 
the same initial eccentricity a to each planet in the system. 
We do not migrate these planets into resonance, but since (in 
general) resonance can add a destabilizing effect ([Chambers | 
let al.|199'6l|Funk et al.|2Ql61|Pu & Wu|2015| ), these simula¬ 
tions can be used as conservative upper limits for the maxi¬ 
mum eccentricity that a system can have and remain stable. 
The results are presented in Figure [bT] and error bars are de¬ 
rived from Poisson statistics. 

Each point in Figure [^represents the fraction of systems 
in our sample having a collision within 2 Myr. We can see that 
for Ci — 0.3, about half of the Kepler systems have had a colli¬ 
sion. Thus from stability arguments, the maximum initial ec¬ 
centricity that planets in our Kepler sample can have in order 
for ~ 50% of them to survive at least 2 Myr is Cmax ~ 0.3. 
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APPENDIX C: PLANET SAMPLE 


Table Cl: Kepler Systems Used In This Analysis 


System 

Planet 

P (days) 

near MMR 

rriplm^ 

rpire 

Np 

M/Mq 

R/Rq 

KOI-142 

b 

10.95 

V 

8.7 zb 2.5 

3.82 zb 0.44 

2 

0.96 zb 0.04 

0.88 zb 0.03 


c 

22.34 

V 

198.8 zb 9.2 

6.82 zb 1.09 




Kepler-120 

b 

6.31 



2.18 zb 0.22 

2 


0.53 zb 0.03 


c 

12.79 

y 


1.53 zb 0.11 




Kepler-127 

b 

14.44 

y 


1.42 zb 0.11 

3 


1.36 zb 0.04 


c 

29.39 

y 


2.62 zb 0.11 





d 

48.63 



2.62 zb 0.11 




Kepler-176 

b 

5.43 



1.42 zb 0.76 

3 


0.89 zb 0.46 


c 

12.76 

y 


2.62 zb 1.31 





d 

25.75 

y 


2.51 zb 1.31 




Kepler-183 

b 

5.69 

y 


2.07 zb 0.87 

2 


0.96 zb 0.41 


c 

11.64 

y 


2.29 zb 0.98 




Kepler-221 

b 

2.8 

y 


1.75 zb 0.22 

4 

0.72 zb 0.05 

0.82 zb 0.07 


c 

5.69 

y 


2.95 zb 0.33 





d 

10.04 



2.73 zb 0.22 





e 

18.37 



2.62 zb 0.22 




Kepler-244 

b 

4.31 



2.73 zb 1.2 

3 


0.8 zb 0.34 


c 

9.77 

y 


2.07 zb 0.87 





d 

20.05 

y 


2.29 zb 0.98 




Kepler-25 

b 

6.24 

y 

9.0 zb 2.4 

2.62 zb 0.0 

3 

1.19 zb 0.06 

1.31 zb 0.02 


c 

12.72 

y 

14.3 zb 2.7 

4.48 zb 0.0 





d 

123.0 


89.9 zb 13.7 

5.46 zb 0.0 




Kepler-267 

b 

3.35 

y 


1.97 zb 0.11 

3 

0.56 zb 0.05 

0.56 zb 0.02 


c 

6.88 

y 


2.07 zb 0.11 





d 

28.46 



2.29 zb 0.11 




Kepler-27 

b 

15.33 

y 

41.8 zb 5.0 

4.04 zb 0.0 

2 

0.65 zb 0.16 

0.59 zb 0.15 


c 

31.33 

y 

21.2 zb 3.2 

4.91 zb 0.0 




Kepler-272 

b 

2.97 

y 


1.42 zb 0.76 

3 

0.79 zb 0.05 

0.93 zb 0.5 


c 

6.06 

y 


1.75 zb 0.98 





d 

10.94 



2.29 zb 1.2 




Kepler-30 

b 

29.33 

y 

11.3 zb 1.4 

3.93 zb 0.22 

3 

0.99 zb 0.08 

0.95 zb 0.12 


c 

60.32 

y 

640.0 zb 50.0 

12.34 zb 0.44 





d 

143.34 


23.1 zb 2.7 

8.84 zb 0.55 




Kepler-305 

b 

5.49 


10.5 zb 2.6 

3.6 zb 0.87 

3 

0.76 zb 0.13 

0.79 zb 0.05 


c 

8.29 

y 

6.0 zb 2.4 

3.28 zb 0.76 





d 

16.74 

y 


2.73 zb 0.44 




Kepler-32 

f 

0.74 



0.76 zb 0.11 

5 

0.54 zb 0.02 

0.53 zb 0.02 


e 

2.9 

y 


1.53 zb 0.11 





b 

5.9 

y 

9.4 zb 3.6 

2.18 zb 0.22 





c 

8.75 


7.7 zb 5.0 

1.97 zb 0.22 





d 

22.78 



2.73 zb 0.11 




Kepler-326 

b 

2.25 

y 


1.53 zb 0.22 

3 

0.98 zb 0.05 

0.8 zb 0.05 


c 

4.58 

y 


1.42 zb 0.11 





d 

6.77 



1.2zb0.11 




Kepler-327 

b 

2.55 

y 


1.09 zb 0.11 

3 

0.55 zb 0.05 

0.49 zb 0.02 


c 

5.21 

y 


0.98 zb 0.11 





d 

13.97 



1.75 zb 0.11 




Kepler-328 

b 

34.92 

y 

28.5 zb 12.9 

2.29 zb 0.98 

2 

1.15 zb 0.22 

1.06 zb 0.44 


c 

71.31 

y 

39.4 zb 13.6 

5.46 zb 2.29 




Kepler-384 

b 

22.6 

y 


1.09 zb 0.33 

2 

0.76 zb 0.05 

0.88 zb 0.25 


c 

45.35 

y 


1.09 zb 0.33 




Kepler-386 

b 

12.31 

y 


1.42 zb 0.76 

2 

0.74 zb 0.05 

0.77 zb 0.43 


© 0000 RAS, MNRAS 000, 000-000 



Tides Alone Cannot Explain Kepler Planets Close to 2:1 MMR 9 



c 

25.19 

y 


1.64 zb 0.87 




Kepler-396 

b 

42.99 

y 

75.5 zb 11.8 

3.49 zb 1.31 

2 

0.85 zb 0.13 

1.06 zb 0.39 


c 

88.5 

y 

17.9 zb 2.8 

5.35 zb 1.97 




Kepler-48 

b 

4.78 

y 

14.3 zb 4.3 

2.18 zb 0.0 

3 

0.88 zb 0.06 

0.89 zb 0.05 


c 

9.67 

y 

9.8 zb 3.3 

3.17 zb 0.0 





d 

42.9 


7.93 zb 4.6 

2.07 zb 0.11 




Kepler-5 6 

b 

10.5 

y 

22.1 zb 3.9 

6.55 zb 0.33 

2 

1.32 zb 0.13 

4.23 zb 0.15 


c 

21.4 

y 

181.0 zb 21.0 

9.83 zb 0.44 




Kepler-5 7 

b 

5.73 

y 

118.1 zb 24.1 

2.18 zb 0.0 

2 

0.83 zb 0.05 

0.73 zb 0.0 


c 

11.61 

y 

7.4 zb 9.4 

1.53 zb 0.0 




Kepler-79 

b 

13.48 

y 


2.62 zb 0.76 

4 

1.1 zb 1.63 

1.4 zb 0.25 


c 

27.4 

y 


2.73 zb 0.87 





d 

52.09 



7.64 zb 1.42 





e 

81.07 



3.38 zb 0.66 




Kepler-81 

b 

5.96 

y 


2.4 zb 0.44 

3 

0.64 zb 0.38 

0.59 zb 0.03 


c 

12.04 

y 


2.4 zb 0.33 





d 

20.84 



1.2 zb 0.33 




Kepler-83 

d 

5.17 



1.97 zb 0.11 

3 

0.66 zb 0.41 

0.59 zb 0.03 


b 

9.77 

y 


2.84 zb 0.44 





c 

20.09 

y 


2.4 zb 0.33 




Kepler-9 

d 

1.59 



1.64 zb 0.22 

3 

1.07 zb 0.05 

1.02 zb 0.05 


b 

19.24 

y 

80.09 zb 4.13 

9.5 zb 0.76 





c 

38.91 

y 

54.35 zb 4.13 

9.28 zb 0.76 
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